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By using spectral methods, we present seven classes of stable and unstable structures that occur 
in a dissipative media. For different parameters of the equation, we find ranges of existence of stable 
structures (spinning elliptic, pulsating, stationary, organized exploding), and unstable structures 
04 (filament, disorganized exploding, creeping). By varying initial conditions, and vorticity, we find a 

reacher behavior of solutions in the form of creeping- vortex (propellers), spinning rings and spinning 
"bean-shape" solitons. Each class differentiates from the other by distinctive features of their shape 
and energy evolution, as well as domain of existence. 
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I. INTRODUCTION 

T— I 

I— I Such as the nonlinear Schrodinger equation (NLSE), the complex cubic-quintic Ginzburg-Landau equation (CCQ- 

^0 GLE) is one of the most intensively studied equation describing weakly nonlinear phenomena in dissipative systems 
Ph [5]- Thus, much work has been done about the features and the properties of this equation as well as its numerous 
^ applications such as nonlinear waves, superconductivity, superfluidity [Bl, Bose-Einstein condensation, Benard con- 
• ^ vection and nonlinear optics 7J. 
^ Solitary waves (or solitons) are self localized solutions of certain nonlinear PDEs describing the evolution of a dissipa- 
tive system [5]. Whereas traditional solitons arc stationary in time and preserve their shape upon interaction, some 
I dissipative soliton solutions of the CCQGLE are non stationary. In Hamiltonian systems, stationary solitons exist as 
^ a result of a balance between diffraction/dispersion and nonlinearity. However, a dissipative system is known to be 
f~ — far from equilibrium and is defined by energy exchanges with external sources. Thus, there is no conserved quantity, 
G\ which implies that they are not Hamiltonian systems. In dissipative systems, solitons are obtained by balancing gain 
Gn and loss. As a result, whereas solutions are defined by one parameter family in Hamiltonian systems, in dissipative 
systems solutions are obtained with varying amplitude and width that are fixed by the parameters of the equation. 
y—^ Therefore, the inverse scattering method [T] which is used to calculate solutions of integrable and some non integrable 
Hamiltonian systems cannot be used in our case. In this paper, we present new classes of soliton solutions of the 
GCQGLE and we analyze both their shape and energy evolution that depends strongly on the parameters used, and 
[ choice of initial conditions. The development of powerful numerical tools over the last decades lead us to use an accu- 
>■ rate spectral method for the spatial discretization, together with a 4*^ order Adams-Bashforth scheme for temporal 
discretization. Due to the intensive computations codes in Fortran 90, we run simulations on ZEUS cluster (256 nodes 
dual Xeon 3.2 GHz processors, f024 KB cache 4GB with Myrinet) at Embry-Riddle Aeronautical University. For 
each class of solitons, simulations lasted from several hours to several days. 



II. NUMERICAL SCHEME OF THE 2D CCQGLE 

Numerical simulations are based on the 2D CCQGLE model ([T]) with transverse Laplacian V^^ that takes into 
account spatial diffraction in the paraxial wave approximation. Quintic terms are essential in order to obtain stable 
solitons [T] 

d,A = eA+ {bl + ici)Viv4 - (63 - ic3)|ApA - (65 - ics)!^!^^. (1) 
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2D solitons 


e 61 ci 63 C3 65 C5 


elliptic spinning and filament 


-0.1 0.1 0.5 -0.8 1 0.04 -0.02 


pulsating and stationary 


-0.045 0.04 0.5 0.37 1 0.05 -0.08 


exploding 


-0.1 0.125 0.5 -1 1 0.1 -0.6 


creeping and propeller 


-0.1 0.101 0.5 -1.3 1 0.3 -0.101 



TABLE I: Initial sets of parameters for 2D solitons. 



A{x, y; z) is the normalized envelope of the field which is a complex function, x and y are the transverse spatial 
coordinates and z is the propagation distance along the optical fiber (or cavity depth). The physical parameters 
of the system are: e linear gain/loss, bi angular spectral filtering, ci — 0.5 second-order diffraction coefficient, 63 
nonlinear gain/loss, C3 = 1 nonlinear dispersion, 65 saturation of the nonlinear gain/loss, and C5 saturation of the 
nonlinear refractive index. The solutions to ([ij will be computed using a Fourier spectral method in which we will 
monitor the energy Q{z) = \A(x,y; z)\'^ dxdy. For a localized solution, Q is finite and changes smoothly 

while the solution stays within the region of existence of the soliton. When Q changes abruptly there is a bifurcation 
and the solution jumps from a branch of solitons that become unstable to another branch of stable solitons, or vice 
versa. As soon as the solution becomes unstable, Q diverges until infinity or collapses to 0. For a certain class of 
solutions, Q will evolve periodically in some regime, and will converge to a finite value [lOj . 

III. PARAMETERS AND INITIAL CONDITIONS 

The main goal of this study is to find continuous parameter ranges around a fixed set in which each class of solitons 
remains stable. 



A. Parameters 

First, an initial set of parameters published in j4] is used for the elliptic, filament, pulsating and exploding classes, 
while for creeping we use parameters published in [5], see Table 1. Starting from these values, for each class we vary 
one parameter in a range of width 0.2 while keeping all the other parameters constant. When we are outside of the 
range of existence of the solitons, we return to the initial value, and vary the other remaining parameters one by one. 
After several simulations, we noticed that 63 and e seemed to be the more sensitive, i.e., small changes of 63 or e 
implied significant changes of the solution. Thus, a step of 0.01 has been chosen for these parameters whereas a step 
of 0.05 has been used for 61, 65 and C5. However, this behavior is not a general rule, so, in some cases, a smaller step 
has been used for the last three parameters. 

B. Initial and boundary conditions 

We use two different types of initial conditions, a Gaussian shape A{x,y; 0) = AqC-'- , and a ring shaped vortex 
A{x,y;0) — Aor"^e~^ e""^ with a vorticity m. Aq is a real amplitude of the initial condition, and 6 = tan~^ (f^) 
is the phase. The inverse widths (shape) of each of the initial conditions are controlled by the parameters ax and 
ay, with r = ^y (axx)'^ + [ayyY, see Fig. 1. Since we used a Fourier spectral method, the boundary conditions are 
periodic and are incorporated in the method. We made sure that the solutions did not reach the boundaries while 
they evolved, otherwise additional reflection would occur, so we increased the domain to avoid reflective waves. 

IV. SIMULATION RESULTS 
A. Elliptic spinning and filament solitons 

1. Description 

For both of these solutions, the same set of parameters is used, see Table 1, line 1. The difference between them is the 
type of initial condition, elliptic or circular. For the elliptic vortex, an elliptic Gaussian pulse with widths ax = 0.15, 
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Parameters 


Initial value Final value 


bi 


0.035 


0.235 


e 


-0.3 


-0.1 


63 


-0.9 


-0.7 


65 


-0.02 


0.18 


C5 


-0.22 


-0.02 



TABLE II: Parameter ranges of spinning stable solitons. 

and cFy — 0.85. This soliton which lacks radial symmetry has two peaks of amphtude diametrically opposed RP. The 
energy first oscillates before converging to a fixed value and remains stationary, see Fig. 2. Increasing the amplitude 
and choosing the initial condition to be circular, Aq = 3.0, Ux = Uy = 0.15 results in a splitting of the initial beam in 
several bell-shaped solitons which are non spinning. As it has been shown in [1] , this class may be unstable and leads 
to chaos, but for some set of parameters, the beams stop splitting and the soliton appears to be stable. This result is 
also observed with higher vorticity. For example, for m — 2 see Fig. 3. 

2. Ranges of parameters and bifurcations 

Elliptic spinning and filament solitons undergo bifurcations which result in stable solitons over a large parameter 
range around the initial sets of parameters of Table 1. By varying one parameter in the range from Table 2, several 
bifurcations happen and other stable solitons appear with an energy evolution similar to these presented in Figs. 
2 and 3. The main bifurcations observed lead to the transformation of the initial beam in a stable spinning ring, 
filament or quite exotic spinning "bean- shape" solitons, see Figs. 4, 5 and 6. 

B. Pulsating solitons 

1. Descnptton 

As it was presented in 0] , a soliton with a very interesting pulsating behavior was discovered for a slightly asymmetric 
{Aq = 5.0, ax = 0.8333, ay = 0.9091) Gaussian beam. For the set of parameters of Table 1, line 2, the soliton evolves 
first as a stationary one, but after a certain time, its energy increases and takes a pulsating behavior before increasing 
again to a new stationary state. During the pulsating period, the soliton keeps a radially symmetric bell shape 
profile at the maximum of energy, whereas for the minimum value of Q, the beam elongates alternatively in different 
directions. However, other pulsating solitons exist and do not evolve in the same way. For example, for two values of 
&3, the energy evolves differently, see Fig. 7. With 63 = 0.37, after a short pulsating behavior, the soliton reaches a 
stationary stable state until z « 420s. Following that a perturbation appears which makes the soliton evolve towards 
a higher state with a new pulsating energy. This pulsating behavior lasts around 100s. For 63 = 0.40, the first 
stationary state does not appear, the energy seems to be collapsing before reaching the pulsating state very fast. 
This pulsating state lasts in this case for almost 200s. Furthermore, the energy values that both solitons reach after 
the pulsating state are relatively different {Q « 60 for 63 = 0.37 and Q w 100 for 63 = 0.40). However, despite the 
differences observed relating to the evolution of the energy, for 63 = 0.40 the structure has the same dynamic behavior 
as described above during the pulsating state, see Fig. 8. 

2. Ranges of parameters and bifurcations 

In the ranges of parameters explored, the pulsating soliton experiences a bifurcation that leads to a stationary bell 
shaped soliton, see Fig. 9. The parameter 65 seems to be a very sensitive parameter for the pulsating soliton. Indeed, 
whereas for 65 e [—0.05; —0.02] the soliton is pulsating, from 0.06 to 0.18 solitons are stationary and from —0.02 to 
0.04 solitons are unstable. Outside of these ranges, from —0.02 until 0.04, solitons are unstable. Table 3 gives an 
overview of the parameter ranges explored for the pulsating and stationary solitons. 
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Parameters 


Pulsating Stationary (non pulsating) 


61 


[0.015; 0.035] [0.085 


0.235] 


e 


[-0.095; -0.045] [-0.245 


-0.105] 


63 


[-0.4; -0.36] [-0.35 


-0.20] 


6b 


[-0.05; -0.02] [0.06 


0.18] 


cs 


[-0.09; -0.08] [-0.23 


-0.11] 



TABLE III: Parameter ranges of existence of pulsating and stationary solitons. 



Parameters 


Range 


61 


[0.135; 0.235] 


e 


[-0.3; -0.1] 


63 


[-1.0; 0.8] 


68 


[0.09; 0.15] 


Cs 


[-0.6; -0.4] 



TABLE IV: Parameter ranges of existence for organized exploding solitons. 

C. Exploding solitons 

1. Description 

For this class we use a radially symmetric Gaussian beam with amplitude Aq = 3, cFx — o'y = 0.3, and parameters 
from line 3 of Table 1. The main feature of this soliton is that it explodes intermittently and periodically resulting in 
significant bursts of energy, while after each explosion, the soliton recovers its initial shape, see Fig. 10. The evolution 
starts from a stable bell shaped structure. Then, circular waves appear around the the initial beam developing from 
center to the exterior. Suddenly and violently, the soliton explodes, it grows up fast keeping radial symmetry, and 
recovers its initial shape after an explosive phase. Even if the soliton always returns to the same shape, the explosions 
are not strictly identical as we can see in Fig. 10. The evolution of the structure during a full explosion process is 
shown in Fig. 11. 

2. Ranges of parameters and bifurcations 

Outside of the ranges of existence of exploding soliton from Table 4, we found two types of solitons. The stable 
stationary solitons, and unstable ones for which the energy diverges to infinity. The later type never recovers their 
shape, and appears probably when the features of the system and the medium (defined by the parameters of the 
CCQGLE) cannot compensate the energy bursts, and hence let the energy grow without bound. 

As it was explained previously, small changes in some parameters dramatically change features of the exploding 
soliton. As it is shown in Fig. 12, the effect of changing 63 shows an energy evolution with random amplitudes and 
shifts. Also, depending on the initial conditions, the explosions may be organized or disorganized, as we can see in 
Fig. 13. 

A ring vortex initial condition has been tested in the same ranges of parameters. Contrary to pulsating solitons 
which did not exist with this initial shape, exploding solitons seemed to be existing. However, even when explosions 
and bursts happened, the behavior of such a soliton was extremely unpredictable. In fact, the evolution of the energy 
is very random, explosions do not appear periodically and the amplitudes of these explosions arc relatively different, 
see Fig. 13. We assume that this initial condition leads to an exploding soliton that is on transition to chaos. 

D. Creeping solitons 

1. Description 

Creeping soliton dynamic behavior also starts from a Gaussian initial condition with radial symmetry, with Aq = 3.0 
and (Jx = (Jy = 0.3. In 2D, creeping solitons behave like tumors. In fact, they spread in all directions until filling the 
domain. Thus, the energy of a 2D creeping soliton is very often characterized by a constant raise, hence precautions 
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Parameters 


Range 


61 


[0.085; 0.135] 


e 


[-0.26; -0.1] 


63 


[-1.0; 0.8] 


65 


[0.13; 0.15] 


C5 


[-0.6; -0.5] 



TABLE V: Parameter ranges of existence for disorganized exploding solitons 



Parameters 


Range creeping Range stationary 


bi 


[0.08; 0.14] not found 


e 


[-0.13; -0.09] [-0.3; -0.14] 


&3 


[-1.3; -1.29] [-1.28 


-1.11] 


65 


[0.3; 0.31] [0.32 


0.36] 


C5 


[-0.101; -0.100] [-0.8 


-0.6] 



TABLE VL Parameter ranges of existence. Left: creeping solitons. Right: stationary solitons. 

must be taken when the size of the domain is defined. Before spreading and fiUing the numerical grid, creeping 
sohtons present a very complex dynamic evolution pattern by taking surprising geometrical shapes that resemble 
fractals. They are still chaotical while but localized in the domain and tend to show some symmetry. In the ID 
case, it has been proved [5] that the region of parameters where creeping solitons exist is filled which a rich variety of 
bifurcations between, stationary, pulsating and creeping solitons. 

2. Ranges of parameters 

Two main bifurcations were observed by varying parameters of Table 1, line 4. The first one takes the shape of 
a ring which is expanding very fast and leads to chaos, the second one leads to stationary solitons. The parameters 
ranges are presented in Table 6. 

By changing the initial condition into a ring vortex of vorticity m = 1, a novel class, spreading vortex soliton 
has been found. This soliton has exactly the same behavior as the creeping, but it is also spinning and that creates 
interesting "propeller" shapes before spreading and filling the domain, see Fig. 15. The energy of such structure starts 
oscillating (during this phase the beam alternatively grows up and comes back to its initial shape) and finally converges 
to a fixed value that takes the shape of a ring vortex. Fig. 16 shows the evolution of the energy of a propeller soliton 
emphasizing this bifurcation: the energy oscillates and seems to start converging but then a perturbation makes the 
energy diverge and increase because of the spreading behavior of the solution. Table 7 shows a comparison between 
the propeller and the ring spinning type. If one compares the energy between the two, it seems like a battle between 
the increasing behavior of the creeping together with the spinning behavior of the ring, when the energy wants to 
remain bounded. 



Conclusion 

In conclusion, we have found fascinating (2-|-l)D structures in dissipative media predicted by the complex cubic- 
quintic Ginzburg-Landau equation. Different types of solutions, spinning or not, such as elliptic vortex, filament. 



Parameters 


Range propeller Range ring spinning 


bi 


[0.08; 0.14] [0.185; 0.235] 


e 


[-0.19; -0.09] [-0.3; 0.2] 


63 


[-1.3 


-1.23] [-1.22; -1.11] 


65 


[0.3 


0.36] not found 


C5 


[-0.1 


-0.06] not found 



TABLE VII: Parameter ranges. Left: spreading spinning solitons. Right: ring spinning solitons. 



FIG. 1: Initial condition. Left: Gaussian shape. Right: ring vortex with vorticity m — 1. 
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FIG. 2: Elliptic spinning soliton. Left: contour plot of the amplitude. Right: energy evolution. 

pulsating, exploding and creeping solitons have been categorized by focusing on their energy, shape and phase evo- 
lution. Regions of existence of these types of soliton have also been calculated in the five dimensional parameter 
space given by the physical parameters of the equation. The study revealed interesting ranges and bifurcations for 
each structure. Whereas non spinning structures have mainly lead to stable stationary solitons, due to bifurcations 
spinning structures, symmetric or asymmetric, split and have revealed to be no spin solitons. By changing the ini- 
tial conditions, and increasing the vorticity, allowed us to observe for the first time new types of solutions such as 
disorganized exploding and propeller solitons. 
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FIG. 3: Filament soliton for m — 2, bs = 0.180 that splits in several bell-shape beams and non spinning phase. Top left: 
contour plot of amplitude. Top right: phase plot. Bottom: energy. 




FIG. 4: Spinning ring soliton with elliptic vortex initial condition m — 1, bg = 
Right: phase plot that shows spinning. 



-0.84, z = 50s. Left: contour plot of amplitude. 
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FIG. 5: Filament soliton with elliptic vortex initial condition m = 1, — 0.18, z — 100s. Left: contour plot of amplitude. 
Right: phase plot that shows no spin. 
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FIG. 6: Spinning "bean-shape" soliton with elliptic vortex initial condition to = 1, e = —0.17, z = 50s. Left: contour plot of 
amplitude. Right: phase plot. 
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FIG. 7: Energy of pulsating solitons. Left; 63 = 0.40. Right; 63 = 0.37. 
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FIG. 8; Shape evolution and elongation during the pulsating state for 63 = 0.40. From top to bottom, z — 90s, z — 120s, z = 
135s, 2 = 150s. 
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FIG. 9: Bifurcations for pulsating soliton for 63 — —0.36 (top energy in red) that becomes stationary bell shaped soliton for 
63 = —0.35 (bottom energy in blue). 
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FIG. 10: Exploding soliton: evolution of the energy which shows periodic bursts and recovery. 
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FIG. 11: Explosion process corresponding to Fig. 10. From top to bottom row wise, z = 10s, z — 12s, z — 14s, z — 15s, 
z = 16s, z = 18s. 
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FIG. 12: Energies of two exploding solitons obtained with 63 — —0.96 (blue) and 63 — —0.81 (red): the explosions are shifted 
with different bursts. 
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FIG. 13: Disorganized exploding soliton starting as a ring vortex with m = 1, and 65 = 0.14 




FIG. 14: 2D Creeping soliton: Gaussian beam is spreading like a fractal. 




FIG. 16: Evolution of the energy of a spreading spinning soliton next to the bifurcation boundary. 



